Làser cooling of trapped Fermi gases 
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The collective Raman cooling of trapped one- and two-component Fermi gases is considered. We 
obtain the quantum master equation that describes the làser cooling in the festina lente regime, for 
which the heating due to photon reabsorption can be neglected. For the two-component case the 
collisional processes are described within the formalism of quantum Boltzmann master equation. 
The inhibition of the spontaneous emission can be overcome by properly adjusting the spontaneous 
Raman rate during the cooling. Our numerical results based in Monte Cario simulations of the 
corresponding rate equations, show that three-dimensional temperatures of the order of 0.08 Tf 
(single-component) and 0.03 Tf (two-component) can be achieved. We investigate the statistical 
properties of the equilibrium distribution of the laser-cooled gas, showing that the number fluctu- 
ations are enhanced compared with the thermal distribution close to the Fermi surface. Finally, 
we analyze the heating related to the background losses, concluding that our laser-cooling scheme 
should maintain the temperature of the gas without significant additional losses. 

PACS numbers: 32.80Pj, 03.75.Ss, 42.50.Vk 



I. INTRODUCTION 

The achievement of Bose-Einstein condensation (BEC) 
p| in trapped dilute atòmic gases has stimulated a large 
interest in the physics of ultracold gases, both bosonic 
and fermionic ones. Related with the latter, similar cool- 
ing methods as those employed to reach BEC have been 
recently employed to accomplish a degenerate Fermi gas, 
i.e. a Fermi gas with temperature T below the Fermi 
temperature (T F ) 0, H S S H • For T < T F the Fermi 
pressure becomes noticeable, and as a consequence the 
Fermi cloud becomes significantly broader than a BEC 
under the same temperature. Additionally, for temper- 
atures below a critical one, T c , the system should un- 
dergo a Bardeen-Cooper-Schrieffer (BCS) transition @. 
When this occurs the system becomes superfluid, due to 
the Cooper pairing of partides near the Fermi surface. 
The accomplishment of BCS pairing and its observation 
have been considered in detail @ ; Efl Ell E3 • However, 
the employed cooling schemes have not yet allowed for 
temperatures below T c , which has been predicted to be 
much smaller than Tp. Those cooling methods rely on 
collisional processes, which in the degenerate regime are 
strongly suppressed due to Pauli blocking. As a con- 
sequence, the cooling efhciency decreases, and technical 
losses have prevented until now from reaching tempera- 
tures well below Tp 0. Recently, however, it has been 
proposed that the latter obstacle could be overcome by 
emp loying Feshbach resonances 0,0] or optical lattices 
|16| in order to increase the ratio T c /Tp up to w 0.2, 
opening the possibility to achieve BCS in a tempera- 
ture range where the collisional processes are still efficient 
enough. In fact, very recently temperatures T < Q.2Tp 
have been reported by O'Hara et al. 0], which are below 
the predicted value of Te- 
la this paper we show that làser cooling could pro- 
vide an effective and realistic tool to cool one- and two- 



component Fermi gases well below Tp. The proposed 
method is realized in the festina lente (FL) regime |l8j . 
in which the spontaneous emission rate 7 is smaller than 
the trap frequeney lu. In this regime the heating intro- 
duced by photon reabsorptions can be suppressed. In 
the absence of technical losses the major limitation to 
the làser cooling in the festina lente regime is set up by 
the finite value of Lamb-Dicke parameter, defined as the 
square root of the ratio of the recoil energy to huj. In 
practice, when the spontaneous emission rate becomes 
much smaller than the trap frequeney, the finite available 
time sets limit on achievable temperatures, due to always 
present loss mechanism. Reaching FL regime is experi- 
mentally feasible, and has been in fact demonstrated in 
experiments by the group of Weiss 0, Híj , in which a 
suppression of reabsorptions for decreasing j/cj has been 
observed for tightly bound àtoms in an optical lattice. 
The làser cooling in the FL regime down to the quantum 
deg eneracy has been already predicted to work for bosons 
plL [2^] and for polarized fermions |2j| . The latter case is 
particularly interesting, since evaporative cooling cannot 
be employed to a single-component fermionic system, due 
to the absence of s-wave scattering,_while less restrictive, 
sympathetic cooling techniques 0, 0, la, > are a l so even- 
tually limited due to Pauli blocking. The BCS transition 
in such a system could be realized by externally inducing 
interatomic interactions, for instance dipole-dipole ones 

In the present paper we analyze the three-dimensional 
làser cooling of Fermi gases down to temperatures T <C 
Tp, by performing Monte Cario simulations of the quan- 
tum dynamics governed by the Master Equation (ME). 
For T comparable to Tp, the làser cooling of fermions is 
affected by the statistical inhibition of the spontaneous 
emission, which results in a decrease of the cooling ef- 
ficieney. In a previous work |2.'"{l| we have analyzed the 
làser cooling of a polarized Fermi gas. We have theoret- 
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ically predicted that the above mentioned problem can 
be overcome by either dynamically adjusting the spon- 
taneous emission rate in a Raman cooling process, or by 
employing appropriately designed anharmonic traps. 

We first consider the single-component case. In partic- 
ular, we analyze the statistics of the laser-cooled gases, 
comparing the results with those expected for a thermal 
distribution. We show that whereas the mean-number 
distribution can be well approximated by a correspond- 
ing thermal one, the fluctuations significantly differ from 
the thermal case. This could have significant effects in 
the BCS pairing of these gases after an intcratomic in- 
teraction is induced in the way discussed above. 

In the second part of this work, we analyze the laser- 
cooling of a two-componcnt Fermi gas. For this case, 
the collisions between different components cannot be ne- 
glected. We consider them within the formalism of Quan- 
tum Boltzmann Master Equation (QBME). This formal- 
ism is shown to produced the expected thermal distribu- 
tion in the absence of làser cooling. We analyze the case 
in which only one of the fermionic species is laser-cooled, 
whereas the other one is sympathetically cooled as a re- 
sult of the collisions between different components. We 
show that the laser-cooling is significantly favored by the 
presence of the second component. We analyzed also the 
statistics of the atom distribution, showing that the fluc- 
tuations at the Fermi surface are significantly enhanced 
compared to the thermal case. Finally, we analyze the 
heating induced by background collisions, and show that 
the làser cooling can be employed to maintain a two- 
component Fermi gas at a fixed temperature for a rela- 
tively long time without significant atom losses. 

The paper is organized as follows. In section II, we 
discuss our cooling scheme, and obtain the ME that de- 
scribes the Raman cooling in the FL regime. We also 
describe in this section the corresponding ME for the two- 
component case. We present the resulting rate equations 
and discuss their physical properties. In section III we 
present our numerical results concerning the làser cool- 
ing of a single component gas, characterizing the statis- 
tical properties of the laser-cooled gas. In section IV we 
present our results for the case of a two-component Fermi 
gas, in which one of the components is laser-cooled. We 
discuss the statistics for this case, as well as the heating 
induced by losses, and its control using laser-cooling. We 
finalize in section V with some conclusions. 



rate 7. The latter can be controlled by modifying the 
coupling from |e) to |r). The àtoms are confincd in a 
non-isotropic dipole trap with freqüències uj^ y ^ z , u x,y,zi 
different for the ground and the excited states, and non- 
commensurable one with another. The latter assumption 
simplifies enormously the dynamics of the spontaneous 
emission processes in the FL limit. When deriving the 
quantum ME, it is important to cancel the off-diagonal 
terms of the density matrix, in order to reduce the ME 
to a rate equation, which is possible to simulate numeri- 
cally with Monte Cario methods. The anisotropy of the 
trap does not need to be very large, and for practical pur- 
poses can be rather small. Its effects in the rate equations 
is not relevant at all, and consequently in our simulation 
we consider the simple case of isotropic trap. The cooling 
process consists of sequences of Raman pulses of appro- 
priate freqüències, adjusted is such a way that they in- 
duce the transition of àtoms to the lower motional states 
of the trap. We assume weak optical excitation, so that 
no significant population in |e) is present. This allows 
to adiabatically eliminate |e), and consequently to con- 
sider only the density matrix p(t) describing all àtoms in 
\g), and being diagonal in the Fock representation corre- 
sponding to the bare trap levels. 

Let us introduce the annihilation and creation opera- 
tors of àtoms in the ground (excited) state and in the trap 
level m (l), which we denote by g m , g^ (e/, ej). These 
operators fulfill Standard fermionic anticommutation re- 
lations: {g m ,gl} = S m ,n, {em^j,} = <5 m ,„. We use the 
Standard theory of quantum-stochastic processes [2^, [2|J 
to develop the quantum ME that governs the dynamics 
of cooling |27j|. The derivation of the ME for fermions 
proceeds basically in a similar way to the corresponding 
derivation for bosons. The latter was presented in de- 
tails in • The ME that describes the dynamics of the 
density matrix p of the àtoms in the FL regime is 



p(t) = Cop + dp, 



(1) 



where 



C p = -iÉ effP {t)+ip{t)Èl ff +Jp(t), (2) 
C lP = -*[£,„, p(t)], (3) 



with 



ei 



II. QUANTUM MASTER EQUATION 

We consider N fermionic àtoms with an accessible elec- 
trònic three-level A scheme, with levels \g), |e) and \r). 
We assume that the ground state \g) is coupled via a 
Raman transition to e) (which is assumed metastable). 
Another làser couples |e) to the upper state |r), which 
rapidly decays into \g). After the adiabatic elimination 
of |r), a two-level system is obtained, with an effective 
Rabi frequeney íi, and an effective spontaneous emission 



i 7^2S,ir,iejg m gl l ei, 



l.m 



(4) 



^ Vim(kL)elg m + H.c, (5) 

(G) 



Jp(t) = 2j^2^ lm gl l eip(t)elg, 

l,m 



Here, 27 is the single-atom effective spontaneous emis- 
sion rate, Í2 is the Rabi frequeney associated with the 
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atom transition and the làser field, uj^ n (uf) are the 
energies of the ground (excited) harmònic trap level 
m (l), 5 is the làser detuning from the atòmic tran- 
sition, £ im = J^d^f^dO s'm8W(9,4>)\i]i m (k)\ 2 , where 
YV(8, (j)) is the fluorescence dipole pattern, and r]i m (k L ) = 

(e,l\e lkL ' r \g,m) are the Frank-Condon factors. In the 
following we assume that Çl < 7 and íi 2 /7 <C uj, which 
allows for adiabatic elimination of excited state |e). To 
assure the FL regime we require also 7 < uj at least in 
the initial phase of the cooling process. 

The dynamics of the two-component system in the FL 
limit, where only the first component is làser cooled, may 
be described by the following ME: 



p{t) = Cop + dp + C 2 p- 



(7) 



The new term £ 2 that appears in Equation JJJ) describes 
the process of elàstic collisions between the two species: 



£-2P 



-i[H co u,p(t)], 



where 



U, 



"■in;q l p9n^mPí9p· 



(8) 



(9) 



Here, b m and b^ are respectively the annihilation and 
creation operators of àtoms of the second component in 
the trap state m. Since at low temperatures collisions 
between àtoms are purely of s-wave type, the operator 
H co ii accounts only for collisions between àtoms in the 
second component and those àtoms in the ground elec- 
trònic state \g) of the first species. The collisions with the 
àtoms in the |e) can be neglected if the làser is sufíiciently 
weak to assure that only a small fraction of àtoms is ex- 
cited in each pulse. The amplitude U n ,m;p,q in Equation 
© reads 



U n , 



4-7T ha* 



m\q,p 



If' 



d 3 xC(x)/4(x)/3 9 (x)</> p (x), 



(10) 

where <f> n ((3 n ) denotes the wavefunction of the n-th level 
of the \g), (|6)) trap, and a sc is the scattering length. The 
term L\ from the ME (JJJ, includes the effective Hamil- 
tonian H e ff 1 which in the presence of two components, 
accounts also for the energies of àtoms of the second com- 
ponent 

Heff = W mSÍSm + ^(^f ~ 5 )4 e l + 

m l 

+ X^» 6 « Ò ™ ~ ^YàrneldmgLei- (H) 

n l,m 

For the assumed regime of parameters, the excited state 
may be adiabatically eliminated. To this end we use 
Standard projection operator techniques j2^. The corre- 
sponding calculus proceed in the same way as in the case 
of bosons, which was described in detail in . Here, we 
present only the final results. For sufíiciently weak in- 
teractions between the two species, the dynamics due to 



the collisions and the dynamics due to the laser-cooling 
arc independent [2é|. The process of làser cooling may 
be described by the ME (JJJ without the collisional part, 
while the collisional part of dynamics after the adiabatic 
elimination takes the form of QBME, similar to the one 
formulated for bosons and presented in |2^,|23|. The rate 
equations for the populations and in level n and 
m of the trap \g) and \b) respectively, are given by 



N 1 



N = 



N 1 ^ 

^n,m p q 



Epopt Ty-l , pcoll 
± p^n ±y p ' / J P-q- 
p m,p,q 

Epopt N i _ pcoii tv 1 ^ 2 (12 

P 

Epcol 
p,q 

n,p,q 

E pcoll 
n.m 



m,p,q 



N 1 ^ 

■*n,m p q 



N^-N 2 

■*p,q n m' 



(13) 



where P°^„ is the rate of laser-induced transition from 
state m to state n of the \g) trap 

D, 2 

popt = /j _ N u 



Ei 



l 2 £,lnVlp( k L) 



[í-(«f-a^)] 2 + 7 a (i2,+6p) : 



(14) 



The rate P^L n , m of a 



with Ri = E„'6n'(l-iVi,) 

collision between two fermions of different species, occu- 
pying the states p and q of the trap respectively, resulting 
in a transition to the states n and to, is given by 



-pcoll 
p,q — *n,m 



m;q,p \ 



^E m ,E p + E q (l-K)(l-N^). 

In the transition rate l|14|l one can clearly identify two 
contributions due to the quantum-statistical character 
of fermions: (a) If the system becomes more degenerate, 
Ri vanishes, and the àtoms remain for very long times 
in the excited state l, i.e. the spontaneous emission is 
inhibited. This affects negatively the cooling process in 
two ways: first, it prolongs the cooling, and second, the 
excited-ground collisions can occur and lead to heating 
and losses. In addition, the adiabatic elimination used 
in the derivation of Equation. (|14f) ceases to be vàlid, 
(b) The fermionic inhibition factor (1 — N n ) appears in 
the numerator of the probabilities, introducing also a 
slowing-down of the cooling process. 

The negative influence of Fermi statistics above can be 
overcome by dynamically changing of 7 in Raman cooling 
|23| . One can increase it gradually during the cooling in 
order to avoid the inhibition effects, but still remaining 
in the FL regime. Note that, the FL condition involves 
the spontaneous emission rate modified in the presence 
of other àtoms: jRi < lo. Still, even if one uses such 
approach some small fraction of the àtoms will remain 
in the excited state after the cooling pulse, and has to 
be removed from the trap in order to avoid non-elastic 
collisions. The latter aim can be achieved by optically 
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pumping the excited àtoms to a third non-trapped level. 
This introduces a new loss mechanism which we take into 
account in the simulations. 

The calculations of the làser cooling were performed 
for the case of a three-dimensional harmònic trap. For 
simplicity we assumed that the trap is isotropic with fre- 
quency lo. The main difhculty for performing Monte- 
Carlo simulations of the laser-cooling for fermions is the 
large number of states in the calculation. Since one state 
can be populated maximally by one fcrmion, in order 
to perform a simulation for a medium-size atòmic sam- 
ple (in our case 10 4 àtoms), we had to consider rather 
large trap ( 80 energy levels in our case). In addition if 
the simulation starts from a temperature above Tp, the 
number of states must be much larger than the number 
of àtoms in the trap. For an initial number of 10 4 àtoms, 
we calculated the cooling dynamics for the trap consist- 
ing of approximately 10 5 states. For a single component 
gas of fermions, where the col·lisions between the àtoms 
are absent, we were able to perform simulations treating 
each state of the trap separately. In this case the tran- 
sition is only induced by the làser, and its rate given by 
H14J depends only on two states. The number of differ- 
ent probabilities that one have to recalculate after the 
change of the atòmic distribution is roughly 10 9 . 

In the case of the two-component Fermi gas, one has to 
account for the process of col·lisions, in which the transi- 
tion rate (|15(l depends on four numbers: two initial states 
and two final states. The number of different probabil- 
ities that one has to recalculate after a single transition 
takes place is roughly 10 18 . The number of operations 
rcquired to calculate the probabilities is orders of mag- 
nitude too large for numerical computations. Therefore 
for simulating the cooling in two-species gas, we have ap- 
plied the ergodic approximation. We assume that the 
populations of the states with the same energy are equal: 
N m = ^E(m)/9E(m), where n G {1,2}, is the num- 
ber of àtoms of the n-th component occupying the energy 
shell E and ge is the degeneracy of the energy shell E. 
This approximation relies on the fact that for typical pa- 
rameters, the collisional processes are much faster than 
the làser cooling. Additionally the thermalization inside 
a single energy shell is much faster than the thermaliza- 
tion between energy shells of different energies. 

The assumption that col·lisions occur much more fre- 
quently than làser induced transitions, is necessary, since 
the làser cooling is not fully "compatible" with the er- 
godic approximation. The cooling is based on the process 
of absorption of a photon from one of three orthogonal 
làser beams, therefore the rate of transition depends on 
the three dimensional structure of the initial and the final 
state. Due to this dependence, we have not been able to 
find a simpler formula, which would be a counterpart of 
Equation l|14|l in ergodic approximation. In the case of 
the collisional process, however, one can further develop 
the rate fTïï|l in order to fully utilize the advantages of 
the ergodic approximation. 

The rate of collision, between fermions occupying the 



energy shells e\ and e2, resulting in a transition to the 
energy shells e% and in the ergodic approximation is 
given by 



P. 



coll 

ei,e 2 — *e3,e4 



ei,e2;E3,e4 



. (16) 

where the amplitude of the transition C/ ei e2 . e3 e4 is de- 
fined as follows 



U, 



ei,e 2 ;e 3 ,e4 



n,m,p,q 
^ | U nra -qp | 



(17) 



In principle, the amplitude of transition U eite2 - e3tei may 
be calculated directly from Equation. (|17J) . by perform- 
ing a summation of the terms | U n>m .q„ | given by formula 
IjlOp. This requires, however, a lot of numerical efforts, 
since the number of operations is of the order N^ t , where 
N st is the number of states in the trap. This can be 
avoided, if the trap is isotropic, in which case we were 
able to derive a simplified exprcssion 



U, 



ei,e 2 ;e 3 ,e4 



2 / A-nha s 



\ to£ 3 



E 



j4 jo 

nmpq nmpq 



(i 2 ) 2 

\ nmpq / 



,(18) 



where £ = y/h/muj. The coefficients Inmpq are defined 
by 



C w = / dxe x x H n {x)H m (x)H p (x)H q {x)i (19) 

J — oo 

where H n (x) denotes the n-th Hcrmite polynomial. Note 
that the summation over indices n (m, p, q) in Equation 
(|18fl . runs from if e„ (e m , e p , e q ) is even, or from 1 
if e n (e m , e p , e q ) is odd, up to e„ (e m , e p , e q ) with an 
increment equal to 2. In order to obtain Equation 118fl . 
we have used Standard summation formulas for Hermite 
polynomials jsj. Equation (|18l) can be further simplified, 
if the collisions involves distant energy shells, i.e. the 
minimum e m j n of the set of numbers {e%, e%\ ez-, 64} is 
much smaller than its maximum e m ax ■ According to Ref . 
[32j, the following approximate formula holds 



U, 



ei,e2;e3,e 4 



/ ATrha s 
\ m£ 3 



(h: 



47T 4 



(20) 



With the help of the exact expression (|18fl , we have ver- 
ified numerically that the approximation 120|) gives cor- 
rect results, already for e max — e m j„ > 2hw. Equation 
(|20|l leads to small errors for low lying energy shells, and 
for the transitions between neighboring shells. In the 
simulations, however, we have used the vàlues calculated 
numerically with the help of the exact formula l|18|) • 



5 



III. RESULTS FOR ONE-COMPONENT GAS 

The numerical simulations of the cooling were per- 
formed using Monte Cario techniques. In the calcu- 
lation we assumed that the àtoms are conhncd in a 
dipole trap, characterized by a Lamb-Dicke parameter 
r\ = 2ira/\ = 2, with a = \Jhj2míü being the size of the 
ground state of the trap, and A the làser wavelength. In 
the context of làser cooling, an optical dipole trap was 
recently employed in experiments on all-optical produc- 
tion of a degenerate gas of two species of Li Q . and it is 
also considered in current experiments in Mg 33] . Due 
to numerical limitations we assumed Lamb-Dicke param- 
eter r\ = 2 [341 ]. The assumed value of rj is not unrealistic, 
and can be realized for instance by trapping potassium 
àtoms in a dipole trap with lu = 2n x 2.4 kHz, which is 
the trapping frequency employed in experiments in Li p| , 
and a làser wavelength A ~ 720 nm. In the calculation 
we use the above mentioned vàlues to calculate real cool- 
ing time. We start our simulations with an initial num- 
ber of àtoms N = 10660, which corresponds to a Fermi 
energy Ep = i'èhw. The three-dimensional isotropic har- 
mònic trap contains 81 energy levels, which corresponds 
to 91881 states. 

In the gas of fermions the main sources of losses are 
background collisions and photoassociation. Background 
collisions result from non ideal vacuum conditions in the 
experiment, and introducc a rclcvant heating mcchanism 
for temperatures well below Fermi temperature j^. In 
the simulations we assumed a background collision rate 
Ibg = 1/350 Hz 01 . Photoassociation losses (when the 
làser is tuned between molecular resonances) are typi- 
cally of the order of 10~ 14 cm 3 /s for làser intensities of 
lmW/cm 2 j3||. In our casc the lascr intensities needed 
are typically 1000 times smaller and we estimate that for 
N = 10660 àtoms, the atòmic density is smaller than 
3.5 x 10 14 cm~ 3 . Therefore the photoassociation losses 
can be safely neglected. The other type of losses which 
we account for in our calculations is the removal process 
of long-living excited àtoms at the end of each cooling 
pulsc. 

Figure n shows the time-dependence of the temper- 
ature of a laser-cooled one-component gas of fermions. 
The inset shows the time dependence of the numbcr 
of fermions, that decreases due to losses. As the ini- 
tial state of the system we consider a thermal dis- 
tribution with temperature To = l.OATp. The cool- 
ing process was divided into three stages, each con- 
sisting of a sequence of two Raman pulses. The em- 
ployed pulses that are used are characterized by the 
following parameters: detuning S/lu — {( — 11, —12), 
(—16,-17), (—19,-20)} respectively, Rabi frequency 
n/7 = {(0.013,0.013), (0.007,0.01),' (0.1,0.8)} respec- 
tively, and length At/u;- 1 = {(250,250), (500,500), 
(2000,4000)} respectively. The use of two pulses with 
slightly different detunings in each cooling stage allows 
to compensate for the oscillatory character of the Frank- 
Condon factors. The vàlues of n were chosen in such way 




[S] 



FIG. 1: Time-dependence of the temperature during the làser 
cooling of a single component Fermi gas by controlling the 
effective spontaneous Raman rate. Inset: time dependence 
of the number of àtoms, which decreases due to background 
collisions and the removal of long-living excited àtoms. The 
different Raman cooling stages discussed in the text extend 
from s to 2 s, from 2 s to 4.6 s, and from 4.6 s to 5.5 s. 




FIG. 2: Time-dependence of the effective single-atom spon- 
taneous emission rate for the case considered in Fig. The 
vàlues for the íirst pulse in the sequence (sòlid line with dia- 
monds) and those for the second pulse (sòlid line with squares) 
are indicated. The different Raman cooling stages discussed 
in the text extend from s to 2 s, from 2 s to 4.6 s, and from 
4.6 s to 5.5 s. 



that not more than 10% of the àtoms is excited during 
each pulse, in order to fulfill the condition of the adia- 
batic elimination. The temperature of the system was 
determined by fitting of the thermal distribution to the 
current distribution of fermions. At the end of cooling, 
the system reaches the temperature T = 0.08Tp. 

The time dependence of the spontaneous emission rate 
is presented in Figure[2 The two curves depict the vàlues 
for the two pulses used at each stage of cooling. From this 
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figure one can see how 7 was controlled in order to speed 
up the cooling, and at the same time avoid inhibition of 
the spontaneous emission, and remain in the festina lente 
regime. In the simulation it was realized by increasing of 
7 gradually during the cooling process, in such a way that 
the average value of the spontaneous emission rate calcu- 
lated for many-atom system: 7^, remains approximately 
constant, and equal to 0.8a;, except for the last pulses. 
In the last cooling stage rather long pulses were applied, 
and therefore 77V could be considered smaller than in the 
other pulses. The expectation value of 7jv is presented on 
Figure|3J It was calculated at the end of each pulse, aver- 
aging over all possible laser-induced transition that may 
take place. With each transition m — ► n one can asso- 
ciate the spontaneous emission rate 7jv(í) = where 
l denote an intermediate state, when the atom is excited. 

At the end of the cooling process, the temperature of 
the gas of fermions stabilizes and the system reaches an 
equilibrium distribution. In the absence of col·lisions, this 
distribution is purely determined by the laser-induced 
transitions, and in general may be very diffcrcnt than 
the thermal distribution. In order to investigate this 
point, we studied the mean population and the fluctu- 
ations of the energy shells. The statistics was calculated 
by time-averaging during the last stage of cooling. To 
this end, the final stage of cooling was extended to the 
time t = 4 x 10 5 üj _1 « 2Qs. This time correspond, in 
average to four absorption and spontaneous emission cy- 
cles for each atom. The average distribution of àtoms is 
plotted in Figure 0] The bars depict the distribution ob- 
tained from the simulation, whereas the sòlid line shows 
the fit of the thermal distribution. According to Fig- 
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FIG. 3: Time-dependence of the expectation value of the 
spontaneous emission rate for the case considered in Fig. 
The spontaneous emission rate 7jv takes into account the in- 
hibition due to the presence of the other fermions in the trap. 
The average was calculated for the atòmic distribution after 
each cooling pulse. The different Raman cooling stages dis- 
cussed in the text extend from s to 2 s, from 2 s to 4.6 s, 
and from 4.6 s to 5.5 s. 




FIG. 4: Mean distribution of fermions (bars) obtained by 
averaging in time during the prolonged last stage of the cool- 
ing process presented in Fig. Q Sòlid line - fit of the thermal 
distribution. 



to 
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FIG. 5: Fluctuations of the distribution of fermions (bars) 
obtained by averaging in time during the prolonged last stage 
of the cooling process presented in Fig. Sòlid line - fluc- 
tuations of the thermal gas, with the temperature and the 
chemical potential determined from the fit of the mean dis- 
tribution. 



ure 0] the mean population numbers are well described 
in the framework of the grand canonical ensemblc. The 
fluctuations of the occupation numbers are presented in 
Figure |SJ The sòlid line as before represents the fluc- 
tuations of the thermal gas, with the temperature and 
the chemical potential determined from the fit of the 
mean distribution. As can be seen, the fluctuations in 
the làser cooled gas have completely different properties 
than the fluctuations in the thermal distribution. In a 
small energy range below Ep, the fluctuations are en- 
hanced, whereas deeply in the Fermi sea they completely 
vanish. From this peculiar behavior we may conclude 
that the làser cooling strongly affects the àtoms near the 
Fermi surface. On the other hand, it does not excite 
those àtoms deeply in the Fermi sea, due to the large 
negative detunings of the employed pulses. 

The existence of such non-thermal fluctuations is a pe- 
culiar feature of the proposed cooling method. The ob- 
served enhancement of the fluctuation near the Fermi 
energy could have an important influencc on the forma- 
tion of Cooper pairs in làser cooled Fermi gases. This 
question, however, lies beyond the scope of this paper 
and it will be the subject of furthcr investigations. It 
is worth stressing that our cooling scheme can be ap- 
plied to a Fermi mixture with scattering length modified 
in the vicinity of a Feshbach resonance. With the help 
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of this technique, the critical temperature for the BCS 
transition is predicted to reach 0.2-0.5 T F 0,0. Làser 
cooling should then allow to achieve such temperatures, 
presumably in a time scale shorter that 1 s (see Figure^l . 



IV. RESULTS FOR TWO-COMPONENT GAS 

In the simulation of làser cooling in the two-component 
Fermi system we considered the same parametcrs as in 
the case of the single component gas. We assumed that 
the two species are placed in the same isotropic optical 
trap with frequency u> — 2ir x 2.4 kHz. The Lamb-Dicke 
parameter is r) = 2 and the wavelength of the cooling 
làser is A ~ 720 nm. The trap considered in the simu- 
lation has the same depth for both species, and contains 
81 energy levels. The losses due to the background col- 
lisions affect in this case both components. In the sim- 
ulations we assumed the same background-collision rate 
Ibg = 1/350 Hz for both species. The two components 
have equal initial number of àtoms N = 10660. The scat- 
tcring length for the interactions between the two species 
is equal to a sc — 157ao, where clq is the Bohr radius. This 
value correspond to the interactions between two spin 
statesLF = 9/2, m F = 9/2) and \F = 9/2, m F = 7/2) of 

4 °k m. 

Our first aim was to investigate the dynamics and the 
properties of the equilibrium distribution of the two- 
component Fermi gas where the col·lisions are described 
by the QBME, and the resulting probabilities l|15fl . In 
this case we do not take into account the làser cooling 
mechanism. The energy levels in the trap are initially 
populated according to a thermal distribution with tem- 
perature T = O.lTp. Then, we let the system evolve in 
the presence of collisions during a time t = 400 ms, corre- 
sponding to about 30 collisions per atom in average. The 
time dependence of the temperatures of the two compo- 
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FIG. 6: Time dependence of the temperatures of the two 
components during the evolution in the presence of collisions 
only. The inset presents the final distribution of the first 
component (upper plot) and of the second component (lower 
plot). 




FIG. 7: Average distribution of àtoms of the first component 
calculated during the evolution in the presence of collisions. 
Sòlid line - fit of the thermal distribution. 




FIG. 8: Fluctuations of the atòmic distribution for the first 
component, calculated during the evolution in the presence 
of collisions. Sòlid line - fluctuations of the thermal gas with 
the temperature and the chemical potential determined from 
the fit of the mean distribution. 



nents is presented in Figure |ü] whereas the inset shows 
their final distributions. As one can notice, the tempera- 
ture in the system slightly fluctuates, but its mean value 
remains constant during the evolution. Therefore the 
final equilibrium distribution can not be very different 
from the thermal distribution. The detailed comparison 
may be carried out on the basis of the mean populations 
and the fluctuations of the energy shells, depicted in Fig- 
ures and [SJrespectively. The bars present the averages 
calculated in the simulation, during the whole time of 
evolution. The sòlid lines depict the corresponding vàl- 
ues for the thermal distribution. We see, that both the 
mean vàlues and fluctuations in the two-component sys- 
tem are very close to the functions calculated for a ther- 
mal distribution. This means that the QBME correctly 
dcscribes the process of thermalization. In addition we 
conclude that for N — 10660 àtoms, finite size effects are 
not observed, and the distributions correspond to the 
ones derived in the thermodynamic limit. 

Let us now discuss the results calculated with the 
inclusion of the cooling mechanism. The time depen- 
dence of the temperature for the cooling of the two- 
component gas of fermions is presented in Figure [ü] 
The inset shows the number of fermions as a function 
of time. The initial state of the system is given by a 
thermal distribution with temperature To = Tp. As 
for the case of single component Fermi gas, the cool- 
ing process was divided into three stages, each consisting 
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of a sequence of two Raman pulses. The pulses that 
are used are characterized by the following parameters: 
detuning S/lü = {(-11,-12), (-16,-17), (-19,-20)} 
respectively, Rabi frequency fl/j = {(0.113,0.113), 
(0.008,0.012), (0.0025,0.004)} respectively, and length 
At/ui- 1 = {(250,250), (2000,2000), (4000,4000)} re- 
spectively. For the considered parameters, not more than 
10% of the àtoms is excited during each pulse. As one can 
observe, a final temperature T ~ 0.03Tp may be reached 
within 6 s. 

We have investigated also the statistical properties of 
the equilibrium distribution in the two-component sys- 
tem, in the presence of both laser-cooling and col·lisions. 
The results are presented in Figure EU (mean occupa- 
tion numbers) and in Figuro 1111 (fhictuations). The 
bars represent the numerical results, calculated by time- 
averaging during the last stage of the cooling. The sòlid 
curves depict the corresponding functions for the ther- 
mal distribution, with the temperature and the chemical 
potential determined from the fit of the average distribu- 
tion. The mean occupation numbers are exactly the same 
as predicted by the thermal distribution. The fhictua- 
tions in the laser-cooled system, however, are substan- 
tially larger than the fluctuations in thermal ensemble. 
The laser-cooling process induces additional migration of 
the àtoms in the region close to the Fermi energy, leading 
to the increased fluctuations in this region. 

Finally we have studied the storage of the degenerate 
Fermi gas produced in the cooling process. The main 
limitation in this case is provided by the background 
losses, which may generate holes deep within the Fermi 
sea. The latter increases the corresponding energy per 
particle, and via collisional thermalization leads to an im- 
portant source of heating [3^ . We calculate the dynamics 
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FIG. 9: Time-dependence of the temperatures of the two 
components (darker and lighter curves) during the cooling of 
a two-component Fermi gas by controlling the effective spon- 
taneous emission rate. Inset: time dependence of the number 
of àtoms. The losses are due to background col·lisions and the 
removal of long-living excited àtoms. The different Raman 
cooling stages discussed in the text extend from s to 1.3 s, 
from 1.3 s to 3.2 s, and from 3.2 s to 5 s. 
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FIG. 10: Mean final distribution of fermions (bars) obtained 
by averaging in time during the last stage of the cooling pro- 
cess presented in Fig. |Ü| Sòlid line - fit of the thermal distri- 
bution. 
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FIG. 11: Fluctuations of the final distribution of fermions 
(bars) obtained by averaging in time during the last stage 
of the cooling process presented in Fig. |Ü| Sòlid line - fluc- 
tuations of the thermal gas, with the temperature and the 
chemical potential determined from the fit of the mean dis- 
tribution. 



of the system, which was aheady cooled to a temperature 
T = 0.03Tp- The time dependence of the temperature is 
presented in Figure^l The plot compares two cases: (i) 
the làser is turned off, and the gas is heated due to the 
background collisions, (ii) the làser is turned on, and the 
cooling pulses of the last stage are continuously applied. 
As we may observe from the figure, in the latter case, the 
làser cooling compensate for the heating induced by the 
creation of holes in the degenerate distribution. Hence, it 
helps to maintain the degenerate gas for a relatively long 
time in the trap. We have also verified that the pres- 
ence of laser-cooling does not lead to substantially larger 
losses. 



V. CONCLUSION 

In this paper we have analyzed the làser cooling of 
trapped one and two-component Fermi gases. We have 
developed the quantum master equation that describes 
the làser cooling of fermions, and the process of colli- 
sions in the case of the two-component gas. We have 
considered the FL regime where the heating due to pho- 
ton reabsorptions can be neglected. We have shown that 
the inhibition of the spontaneous emission due to Pauli 
blocking can be overcome employing collective Raman 
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FIG. 12: Heating due to background collisions in the degener- 
ate two-component Fermi system. The background collisions 
rate is assumed to be 7 = 350 Hz. The temperature of the 
system slowly increases when the cooling làsers are turned off. 
In the presence of the cooling làsers the temperature of the 
system remain unchanged. 



cooling with dynamically modified spontaneous Raman 
rate. Additionally, a slightly anharmonic trap could be 
cmployed. By means of Monte Cario simulations of the 
corresponding rate equations, we have shown that it is 
possible to cool one- and two-component Fermi gases 
down to temperatures of the order of a few percent of 
Tp. The results obtained in this paper present by no 
means fundamental limits of làser cooling; temperatures 
below 1% of Tp can be achieved by optimizing the cooling 
procedure |38| . 

The làser cooling is an external process that influences 



the dynamics of trapped àtoms, therefore one can expect 
that the equilibrium distribution of the laser-cooled sys- 
tem may differ from the thermal distribution. We have 
analyzed, both for the single- and the two-component 
case, the statistics of the laser-cooled gas, comparing the 
results to the predictions of the grand canonical ensem- 
ble. We have shown that although the mean-number dis- 
tribution can be well approximated by a thermal one, the 
fluctuations are significantly different. In particular, the 
fluctuations are considerably enhanced at the Fermi sur- 
face. This effect could play a significant role in the BCS 
pairing in làser cooled gases, and it will be the subject of 
further analysis. 

Additionally, we have discussed the possible loss sour- 
ces that may affect the làser cooled gas. In particular, 
for the two-component gas, the creation of holcs deep in 
the Fermi sea due to background collisions could lead to 
a significant heating of the sample, which should impor- 
tantly reduce the life-time of the degenerate gas. In this 
sense, we have shown that our laser-cooling scheme can 
be employed to maintain a two-component Fermi gas at 
a fixed temperature without significant additional atom 
losses. 
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